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Abstract. We establish the global existence of L°° solutions for a model of 
polytropic gas flow with diffusive entropy. The result is obtained by showing 
the convergence of a class of finite difference schemes, which includes the Lax— 
Fricdrichs and Godunov schemes. Such convergence is achieved by proving the 
estimates required for the application of the compensated compactness theory. 



1. Introduction 

We consider the following system modeling isentropic gas flow with smoothly 
varying entropy. The model reads in Eulerian coordinates 

(1.1) Pt + m x = 0, 

9 

772 

(1.2) mt + (— +P(P,S)) X = 0, 

P 

(1-3) (pS) t + (mS) x = (~S X ) X , 



P 



where 



p(p,S) = Ke^ s / m p^, 

where *H > and 7 > 1 are constants, and n = ^(7 — l) 2 . Here p represents the 
gas density, m is the momentum defined as m = pu, where u is the gas velocity, 



p represents the gas pressure, and S stands for the entropy. The system (1.1 1 



(1.3) is a mathematical model intended to approximate the more physical model 



where equation (1.3 1 is replaced by the energy conservation law, which for smooth 
solutions is equivalent to the equation (pS) t + (mS) x = 0, and this motivates our 
mathematical model. 

Initial data are given by 

(1.4) p{x,0) = po(x), m(x,0) = m (x), 

(1.5) S(x, 0) = S (x) = a(y (x)), y (x) = / p (z)dz. 

Jo 

Assume that 

(1.6) Po, mo,— eL°°(R), Po >0, <reW^(R). 

Po 
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In particular, the initial data (and the solution) allows for the occurrence of vacuum. 
In addition, we also assume that a is periodic with period, say, 2ir, that is, 

(1.7) a(y + 2n)=a(y), y G K. 



We remark that assumption (|l.6|), imposed on a, implies that the solution of the 

-(y-^ 2 /^a(z)dz, 



heat equation with initial data a, 



(1.8) 
satisfies 



a(y,t) := 



1 



(47rt)V2 



(1.9) \&(y,t) - <j\,\a y {y,t)\,\&yy(y,t)\ < C e \ with a :- 



1 

2^ 



2 71 



a(z) dz, 



for some absolute constant Co > 0. Indeed, (1.6) and (1.7) imply the absolute con- 
vergence of the Fourier series of a, a' and a". On the other hand, a straightforward 
calculation shows that 

1 r , ^2,, ., (-y 2 + (y+2iktf)/U 

I e -(y- z f/4t e ikz dz = e 



(Ant) 1 / 2 



(4 7 rt) 1 /2 

ihy — k 2 t 



e 4t dz 



for any fc G K, which then gives the asserted asymptotic behavior, by plugging the 
Fourier series for a, a' and a" in (1.8) and the corresponding equations for a y and 
a yy , obtained from (1.8) by replacing a by a' and a", respectively. 
We have the following definition of weak solution. 

Definition 1.1. We say that a function (p,m,S) £ L°°(M. x (0,oo)) is a weak 
solution to p3|)-([L5|) if: 

(i) m/p G L°°(R x (0,oo)); 

(ii) for all 4> G C£°(R 2 ), 



(1.10) 

(iii) 
(1.11) 

(1.12) 



{p,m)(x,t)4>t + (m, \-p(S,p))(x,t)4» x dxdt 

x(0.oo) P 



+ / ( / o ,m )(x)^(x, 0)dx = 0; 



where 



S(x,i) = 



y(x,t) 



/47Tt 



e 4t cr(2) ate, 



p(z, f) dz 



m(0, s) ds. 



We observe that away from vacuum, equation (1.3), through the Lagrange trans- 
formation, (x,i) i y (y(x,t),t), with y(x,t) given by (1.12), becomes 



St — Syy, 

and this justifies (iii) of Definition |1.1| 

Indeed, the interplay between the Lagrangian and Eulerian formulation of the 
model is important. For the record we note that the model ( 1.1 )— ( 1.3 1 reads in 
Lagrangian coordinates 

V t ~Uy= 0, 



(1.13) 



U t + P(V, S)y 



0. 



Sf — S; 



<)<) 
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where v = 1/p is the specific volume. We remark that, despite the fact that system 
(1.13| has a form much simpler than ( |1.1[ ),( [L2| ,( 1.3 1, the possibility of occurrence 
of vacuum turns the direct analysis of the Cauchy problem for ( 1.13 ) a very difficult 
task and so, as in the isentropic better strategy is to proceed with the the 

analysis of the corresponding problem in Eulerian coordinates, that is, ( |l.l| -( |l~5] ). 
Our main result reads as follows. 

Theorem 1.1. There exists a constant r( 7 ) > such that if \\(po, m o)l|oo < r {l)j 
then there exists a global weak solution to the Cauchy problem ( 1.1 )- |L5| ) satisfying 
an entropy inequality of the form 

(1.14) r) :t (p,m,S)t + q*(p,m,S) x < -Ce~*, 

in the sense of distributions, for some C > depending on L°° bounds for p, to, S, 

where 

(1-15) 

?7*(p, to, S) = -pu 2 H eyi-VS/M-pi ) q*(p,m,S) = un*(p,m,S) + pu. 

i 7 — 1 

Moreover, r( 7 ) — ¥ oo as 7 — > 1+. Further, if po,mo are periodic with period L such 
that yo(L) = 2n, we have the following decay 

r L 

(1.16) lim / \(p(x, t), m(x, t), S(x, i)) — (p, to, S) \dx = 0, 

t-H3Oj 

where p,rh,S are the mean values of po, Too, So, respectively. 

2. Background results 

Let us first recall results for the p-system for a polytropic gas in Eulerian coor- 
dinates. More precisely, we consider the system 

(2.1) p t + m x = 0, 

9 

777, 

(2-2) m t + (—+p(p)) x =0, 

where the pressure is given by p(p) — Ke^ 7 ~ 1 - )S '/ 5R p 1 . For later use we observe that 
we can rewrite the conserved quantities in terms of the other variables, viz., 
(2-3) 

p = p{p, S) = (^) 1/7 e (7-DS/(7«) ; m = m{U}P: S ) = pu = u(V) lh e^ s l^\ 

Here we consider the isentropic case where the entropy S is considered a constant. 
Recall that the functions 



(2.4) W = u+U Pp y / =u + e es /*p< 



(2.5) z = u -±( Pp y« = u - e es/* p e = u _ { y ) 



with = i(7 — 1), form a pair of Ricmann invariants for system (2.1 )-( |2T2| in the 
isentropic case where S is constant. A standard calculation (see, e.g., [61 [3]) yields 
that the rarefaction curves are given by 

to = %± 7 vv^(/-A 
Pi 



while the Hugoniot locus reads 

to = % ± ee° s i*p(±{pi - p]){p - Pl )) 1/2 , 
Pi Wi > 
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2.0 2.5 



Figure 1. The wave curves for the p-system for a given left state. 



from a given left state (pi,mi). When we involve the entropy condition we find that 
the wave curves equal 

(2.6) 



mi 
m = — p 

Pi 



a es/<R 



mi 
m = — p 

Pi 



(2.7) 

W 2 (pi,mt) 

In the variables (p, u) we find 

(2.8) Wi{p u ui): u = m- 

(2.9) W 2 (pi,ui): u = ui + 



P^iP" - Pl)iP - Pi) 



o^p^y-pDip-pi) 
' y/v s /°v - Pi) 

ee^[^-pJ)(p- Pl ) 
ee^{^-pj){p-pi)) 

7 l/2 e 0S/!R( p _ ^ 



1/2 



1/2 



1/2 



1/2 



for p < pi, 
for p > pi, 

for p < pi, 
for p > pi. 

for p < pi, 
for p > pi, 

for p < pi, 
for p> pi- 



An important property of the p-system is that the Riemann invariants provide 
invariant regions. More specifically, (see, e.g., [3J Lemma 5]) if (p (x), m (x)) € 
SI = {(p,m) | w < wo, z > z , w — z > 0} for all x € K, then also the solution 
(p(x, t),m(x, t)) will remain in SI, that is, (p(x, t), m(x, t)) £ SI for (x, t) € fx [0, oo). 

An entropy-entropy flux pair (77, q) for the p-system satisfies for smooth solutions 

r)(p,m) t + q{p,m) x = 0. 



Consistency with the system (2.1)- (2.2) requires 

(2.10) Vq(p, m) = Vr?(p, m)VF(p, m), 

where F — (m, ^- +p(p)) is the flux function of the p-system. A particular choice 
of entropy-entropy flux pair (77* , q* ) reads 



(2.11) 
(2.12) 



2p 7 -l 



-pu 



7-1 



= ( 7 -l)S7-R n 7 



q* = urj* + pu — urj* + une^ 1 1 ^ s ^ R p 1 . 



More generally, the weak entropy-entropy flux pairs (77, q) constitute a class of 
entropy-entropy flux pairs of particular interest in isentropic gas dynamics, as first 
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Figure 2. The triangle in the Riemann invariants (left) is mapped 
into the indicated region bounded by rarefaction waves in the 
(p, m)-plane. 




Figure 3. The solution of the Riemann problem using the Rie- 
mann invariants as coordinates. Through the right state backward 
Riemann invariants are drawn. 



pointed out in [5], and they are characterized by the following conditions at the 
vacuum line: 

r](p,u)\ p=0 = 0, i] p (p,u)\ p=0 = g(u), 
for some continuous function g. Let us denote 

X (p,u;S) = (P-u 2 ) x +1 
P 

where (a) + = max{0,a} and A = " ^ s observed in [8], weak entropy-entropy 

flux pairs can be given by the integral formulas 

(2.13) »?(P,«)= / ff(0x(ft£-«K. 



(2.14) 



q(p,u)= / g(0(6Z + (l-6)u) X (p,t;-u)dt. 



Remark 2.1. Observe that the entropy pair (rj*,q*), defined in (1.15), is a weak 
convex entropy pair. Moreover, for any weak entropy pair (77, q) there exists a 
constant C v > such that rj + C^ry* is convex. 
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Let us now turn to the full system 
(2.15) p t + m x = 0, 



m t + {— + p(p,S)) x = 0, 
P 



(2.16) 

(2.17) {pS) t + (mS) x = 0, 

where the pressure p is given as above. The Riemann problem is the initial value 
problem for the system ( 2.15 )-( 2.17) with special initial data consisting of a single 



jump between two constant states, viz. 



(x) 




for x < 0, 



for x > 0. 



The system ( 2.15 )-( 2.17) possesses three eigenfields associated with the eigenvalues 
Xi=u - y/p~ p , A 2 = u, A 3 = u 



Pp- 



The solution to a Riemann problem for system (2.15 )-( 2.17 ) may be described using 
the coordinates w, z, S, that is, the Riemann invariants for the p-system and the 
entropy, in the following way. Consider first the case when the solution does not 
contain vacuum. The solution of the Riemann problem, starting from the left state 
(pi, mi, Si), consists of a slow wave in which the entropy S remains constant (i.e., in 
the (w, z)-plane determined by S = Si), followed by a contact discontinuity in which 
the velocity u and the pressure p remain unchanged, and finally a fast wave with 
constant entropy S (i.e., in the (w, z)-plane determined by S = S r ) connected with 
the given right state {p r ,m r , S r ). Along the slow wave we can write the Riemann 
invariants a|3 

w = u 1 {p;pi,ui,Si) + e eSl ^p e , 
z = u 1 (p;p l ,u l ,S l )-e 6S ^p e , 



where u = u\{p', pi,Ui,S{) is the slow wave given by (2.6). For the fast wave we 
consider the backward wave (i.e., consisting of the states that can be connected to 
a given right state from the left), and the Riemann invariants read 

w = u 2 (p; Pr,u r , S r ) + e 9Sr / m p 6 , 

z = u 2 {p; p r ,u r , S r ) - e eSr/ ^p 6 , 



(2.19) 



where u = U2(p; p r ,u r , S r ) is the fast backward wave corresponding to (2.7). The 
contact discontinuity, with pressure p* and velocity u* , jumps from a left density 
Pi to a right density p* determined by 



(2.20) 

which yields 
(2.21) 



P 

* 

u 



(ft* 



,(-y-l)S r /SH 



(PrY' 



Ul(P*]PhUl,Sl) = U2(p*;Pr,U r ,S r ), 



P± = e (7-l)(Sr 
p* 



-Si)/(7iH) 



4t turns out to be easier to describe the solution using the speed u rather than the momentum 
m as a variable. 
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z 




Figure 4. The slow Riemann invariant through the left state (blue 
curve) , and the backward fast Riemann invariant through the right 
state (red curve). In addition the yellow curve (w,z), whose in- 
tersection with the slow Riemann invariant determines the contact 
discontinuity. 



to be inserted in the second equation for the velocity, u\ =112, to determine p\ 
and p*. In terms of the Riemann invariants we find that w jumps from u* + 
{p* / K ) e h e - es '/^ to u* + (p* /n) e ^e- es -/ { ^\ and similarly z jumps from u* ~ 
(p* /K) e h e - es il^) to u* - (p*/K) ^e- eSr / < -' f ^l An alternative way to describe 
the contact discontinuity is the following. Consider a point on the backward fast 



wave curve with Riemann invariants {w,z) given by (2.19), which we can write as 
w = u 2 + (p/K) e ^e~ es -/^ and z = u 2 - {p/ K) e r^ s -/^ . Construct now 
another curve (w, z), given as a Riemann invariant with the same velocity U2 and 
pressure p as (id, z), but with the entropy S r replaced by Si, that is, 

w = u 2 + (p/ K ) 6 /-'e- es '/^\ z = u 2 - (p/nf^e-^/W. 

We find 

W + Z = 2ll2 = w + z, 

w-z = 2(p/ K ) e ^e- es ^^ =(w- z)e e ^- s ^/^\ 



which yields 



, r = H(l + e 0(S,-S r )/(7«)) + £(1 _ J&SrVWl)^ 



2 v ' 2 

Z=E(l- e e(Si-S r )/( 7 !H)) + |(1 + e 0(S,-S r )/(~f<R)y 

The intersection between the slow wave curve in the Riemann invariants plane and 
the curve (w, z) determines the values of the variables to the left of the contact 
discontinuity. Through this intersection we draw the line where w + z is constant, 
and the intersection between this line and the backward fast wave gives the values 
of the variables to the right of the contact discontinuity, cf. Figures [4] and [5] 
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Figure 5. The same data is in Figure [4] Curves for the invari- 
ant region for the corresponding p-system are added (black). In 
addition, the dashed line is given by w — z equals a constant de- 
termined by the intersection between the yellow and blue curves. 
The interaction of this straight line with the red curve gives the 
value on the right of the contact discontinuity. The right figure is 
a close-up near the intersection. 



The solution involves vacuum when the slow wave is a rarefaction wave that 
connects to a state on the vacuum line w — z; the velocity is then given by u* = 
Ul + '-y 1 / 2 e e ' s '/ K pf and w — z = u* . Similarly, the given right state connects 
via a rarefaction from a vacuum state with velocity u* = u r — r y 1 ^ 2 e eSr ^pf and 
w — z = u* . 

3. Proof of Theorem 11.11 

3.1. Construction of approximate solutions. Here we provide the full proof 
of Theorem |1.1| We construct approximate solutions for ( 1.1 )— ( 1.3 ) by using a 



Godunov-type finite difference scheme based on solving Riemann problems at each 
time step, updating the approximate S using the Lagrange transformation, and 
averaging at the end of each time step. 

Before we begin the proof, let us describe the fundamentals of the construction 
of the approximate solution. We discretize both in space and time. Let h = At, 
and Ax = ch with c > to be chosen by the CFL condition 

c > sup m h f X '*} ± Jp P {p h (x,t),S h {x,t)) , 

(x,t)e]Rx[0,oo) P \ x 7 l ) v 

which is possible as long as we can obtain an L°° a priori bound for 

The initial data po,mo,So is approximated by step functions with jumps at 
x i-i/2 '■= (* — l/2)Ax for i 6 Z. The multiple Riemann problems are solved for 
t e [0, h). At t = h a new step function is created with jumps at Xj-i/a (details 
given below), and new Riemann problems are solved. More precisely, suppose the 
approximate solution U h = (p h ,m h ,p h S h ) has been defined for t < jh and that 
U h (x,jh) is constant for x € I, where 



Ii = (x i _ 1 / 2 ,x 4+ i/ 2 ), iez. 
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For t £ [jh,(j + l)h), setting Xi — iAx, i £ Z, we define U h {x,t) by glueing 
together the solutions of the Riemann problems for the system (2.15 1— (2.171 de- 



fined at [xi, Xj+i] x [jh, (j + l)h), determined by the discontinuities at the points 
{ x i+i/2ijh), i G Z. Inductively this yields a function U h defined on M. x [0, oo), as 
long as we are able to obtain the necessary a priori bound mentioned above. 

We now provide the details of the construction of the approximate solution. 
Assume that we have constructed the approximate solution U h for x £ R and 
t < jh, and have defined it at time t = jh as a piecewise constant function with 
jumps at x i+ i/2 for i G Z. For (x,t) £ [xi,Xi + {\ x [jh, (j + l)h), i £ Z, let U h (x,t) 
be the solution of the Riemann problem (2.15 1-( 2.17) as described in the previous 
section. Set 

y h {x,t)= ( p h {z,t)dz- ( m h {0,s)ds, x G K, t G [jh, (j + l)h), 
Jo Jo 



and 

a h (x,t) - 

We then defin^3 
(3.1) 

(3.2) 

(3.3) 

for x £ Ii. 



' ' e~ iv (a: « *' a{z)dz = a(y n (x,t),t), (x,t) G K x [0. (j + l)h). 
Ant ' 



p h (x, (j + l)h) = / p h (x, (j + l)h - 0) dx, 
Ax Jh 

m h (x, (j + l)h) = ^f ™' l (i, (j + l)h - 0) dx, 



S h (x, (j + l)h) = -i- / a h (x, (j + l)h - 0) d£, 



Ax 



3.2. Convergence proof. We now address the questions of the L°° a priori bound 
and convergence of U h as h — > 0. First, we investigate the problem of obtaining an 
a priori L°° bound for the approximate solution. Let us denote 

w h (x, t) = w(U h (x, t)), z h (x, t) = z(U h (x, t)). 

Let r > be such that 

w h (x, 0) < r, z h (x, 0) >-r, x £ R. 

We assume for the moment that w h , z h satisfies an a priori bound of the form 

(3.4) w h (x, t) < R, z h (x, t) > -R, (x, t) £ K x [0, oo), 

for some constants R > r, and we will find a condition relating r and R under 
which (3.4) can be justified. 

We first observe that if ( |3.4[ ) holds, then, for any (x\,ti), (2:2,^2) £tx [0, 00), 

(3.5) \y h ( Xl , h) - y h (x 2 ,t 2 )\ < C(R)(\ Xl - x 2 \ + \h - t 2 \ + h), 

for some constant C(R) > depending only on R. In what follows, C(R) will always 
represent a positive constant depending on R that may differ from one occurrence 
to the next one. 



2 We use the standard notation f(x ± 0) = lim e ^o f( x i e). 
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Figure 6. Assuming that the initial data are in the shaded region, 
we show the existence of an R such that the solution remains in 
the larger triangle. The vacuum line is w = z. 



We also observe that 

\a h (x u t 1 )-a h (x2,t 2 )\ 

= \d{y h (x 1 ,t 1 )M) ~ a(y h {x 2 ,t 2 ),t 2 )\ 
< \a{y h {x l ,t l ),t l )-a{y h {x2,t2),t x )\ 

+ \a(y h (x 2 ,t 2 ),t 1 ) - a(y h (x 2 ,t 2 ),t 2 )\ 



(3.6) < \y h (x x , h)-y h (x 2 ,t2)\ / \u y (ry\ + (1 - r)^,^)] dr 

Jo 

+ |ti-*a| / \a t (y^,et 2 + (1-6)^)1(19 



< C(R) - x 2 \ + |ti - t 2 \ + h)e-^ + \h - t 2 \e-™< t ^) 

< C(R){\xi -x 2 \ + \h -t 2 \ + h)e- min ^> t2 \ 



where we have used (1.9) and denoted = y h (xi,t i ) 1 i = 1, 2. 
Assume inductively that 

w h {x,t) < rj, z h (x,t) > ~rj, (x,t) eRx [0,jh], 

for some constant Tj. For t 6 [jh, (j + l)h) the approximate solution is defined by 
solving the Riemann problems given by the discontinuities at the points {xi + i/ 2 , jh), 
i £ Z. Since the p-system enjoys an invariant region given in terms of w and z, the 
only possible increase in w beyond rj , and, similarly, the only possible decrease in z 
beyond — rj, may occur across the contact discontinuity. Here both the velocity and 
the pressure remain unchanged, and the sole change is in the entropy. Observe first 
that since the slow Riemann invariant is increasing in w, there can be no increase 
in the value of w. Fix x and let t £ [jh, (j + l)h). We see from Figure [7] that 
the vertical line x equals a constant crosses slow or fast waves before it crosses the 
contact discontinuity. Let jh < t < t < (j + l)h denote two times such that t is 
after the fast or slow wave, but prior to the contact discontinuity, while t is after 
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(i-3/2) Ax 



(i- 1/2) Ax 



(i+1/2) Ax 



Figure 7. Schematic figure of the solution of the Riemann prob- 
lem. Contact discontinuities are indicated by thick lines. We see 
that the vertical line x equals a constant first intersects a slow or 
a fast wave before it crosses the contact discontinuity. 



the contact discontinuity. Then we find 

z h {x,t) = z h (x,t) + (z h (x,t)- z h {x 1 i)) 

> z h (x,t) - \z h (x,t) - z h (x,t)\ 

> -Tj - \z h (x,t) - z h (x,t)\ , 

since the solution of the p-system remains within the invariant region. Furthermore, 
\z h (x,t) - z h {x,t)\ 



< 



(P^-sV*)^ t) - (?-e- sh ^) 6 ^(x, t) 

K K 



< (?- e - sk ^) e Mx, t) e es "W> r t) 

K 



(3.7) 



x 



r 9sh /^(x,i) - e- esh ^(x,t) 



< l -( w h ~z h ){x,t)e esh '^Hx,t) 



-es h /{~tm) 



(x, t) — e 



-es h /( 1 m) 



(x,t) 



< rj C(R)( 7 -l)\lS h (t)}\ 
= rj C(R)( 7 -l)\lS h {jh)j\ 



as we have replaced both the jump in the exponential by the corresponding jump 
in the exponent and estimated e es /^^{x,^) by a common constant C(R). Next 
we estimate the jump in the entropy. Let x\ and X2 be two points on the left and 
right side of a jump, respectively, thus x\ < Xi-1/2 < £2, with X2 — x\ < Ax. We 
obtain 

\[S h (jh)]\ = \S h (x 2 ,jh)-S h ( Xl ,jh)\ 

' \a h (x + Ax,jh) - a h (x,jh)\ dx 



(3.8) 



< 



< C{R)he 



-jh 



by (j3_6|. 

This yields 

(3.9) z h (x,t) > -r i (l + C(i?)(7-l)/ie-^), 
and we conclude that 

(3.10) z h (x,t)>-r J {l + C{R){ 1 -l)he-i h ), t e \jh, (j + l)h). 
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A similar calculation leads to 

(3.11) w h (x,t)<r j (l + C(R)( rf -i)he-1 h ), t € \jh, (j + l)h). 



At i = (j + l)h we average the approximate solution as described in (3.1 1 



(3.3). Here we argue as follows. We first observe that the averaging of the values of 
(p h (x, (j+l)h-0),m h (x, (i+l)/i-0)) in the intervals := h x {t = (j+l)h-0}, 
i £ Z, in order to obtain the values of (p h (x, (j + l)h),m h (x, (j + l)/i)) in these 



intervals, does not affect the bounds (3.10) and (3.11 ). More precisely, at each such 



interval, S h (x, (j + l)h — 0) assumes at most 3 values, due to the possibility that two 
contact discontinuities, departing from (x^i/ 2 , jh) and {x i+1 / 2 , jh), respectively, 
end inside Z| + . This means that the values of (p h , m h ) in each interval I? + belong 
to the union of at most 3 regions of the form 

R a := {(p, m) : -Gp + e es ^p 6+1 <m<Cp- e es ^ m p e+1 }, a =1,2, 3, 

for some constant C > common to all regions R a , a = 1,2,3. But, one easily 
check that Si < S 2 implies R\ D R 2 , that is, the regions R a , a — 1,2,3, are 
contained in that one corresponding to S 1 * = minj^i, S 2 , S3}. In particular, if we 
define 

S*(x, (j + l)h) := mm{S h (Z, (j + l)h - 0) : £ e h}, in i e Z, 

then, from the convexity of the regions R a , we have 

(3.12) z(p h ,u h ,SZ)(x,(j + l)h) > - rj (l + C(R)^-l)he^ h ), 
and also 

(3.13) w(p h ,u h ,S!:)(x,(j + l)h)<r ] (l + C(R)( 1 -l)he^ h ), 



where 



I u h (x, (j + l)h — 0), otherwise 

and we agree that the value of u h (x, (j + l)h — 0) at a vacuum interval is the mean 
value between its values at the extremes of the interval, which determines precisely 
the values of u h (x, (j + l)h — 0) for all i£i Observe also that the case in which 
Ij + is contained in a vacuum interval is trivial since p h = m h = in such an 
interval, and so the values of p h and m h do not change through averaging on . 



Now, we need to check how the bounds ( 3.12 ) and ( 3.13 ) change when we replace 



S*{x, (j + l)h) by the values of S h (x, (j + l)ti) given by (3.3). For this, we first 



estimate the change in S h from S h (x, (j + l)h-0), to S h (x, (j + l)h), given by p3| 
As already mentioned, S {x, (j + l)h— 0) can be one of three values; either the value 
S h (x,jh), or the values of S in the neighboring intervals, that is, S h (x ± Ax,jh). 



In any of the three cases, the entropy is given by a formula similar to (3.3), but 
with (j + l)h replaced by jh. We consider the most representative case where the 
value is in a neighboring interval. Thus 

\S h (x,(j + l)h)-S h (x-Ax,jh)\ 
(3.14) -AT I \<r h &U + 1 ) h )-° h (%- &x,jh)\ dx 



< C(R)he 



again by (3.6). Since 



z h (x, (j+l)h) = z(p\ u\ S h )(x, (j+l)h), w\x, (j+l)h) = z(p h ,u h , S h )(x, (j+l)h), 
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we conclude as above that 

z h (x,t)>z(p h ,u h ,S*)(x,(j + l)h) 

- \z{p\ u h , S*)(x, (j + l)h) - z(p\ u\ S h )(x, (j + l)h)\ 
>-r i (l + Cr(J2)( 7 -l)^e-^) a =:-r i+lj 
w h (x,t)<w(p\u h ,S>?)(x,(j+l)h) 

+ \w(p h , u\ S*)(x, (j + l)h) - w{ P \u\ S h )(x, (j + l)h)\ 
< rj (l + C(R)(j-l)he-i h ) 2 =r j+1 . 



(3.15) 



It remains to estimate the r,. From the inductive formula (3.15) for the r 7 *, we 



find 



-kh\2 



(3.16) 



= r]l(l + C(R)( 7 -l)he- 

fc=l 

j 

< rexp (2C(i?)( 7 - 1) ^ e - kh h)) 



k=l 



< rexp(2C(i?)(7 

< re 2C(K)( 7 -l)_ 



s ds) 



Therefore, we see from (3.16) that the condition relating r and R under which 



the a priori bound (3.4) holds is 
(3.17) 



Jte-2(7-l)C(fl) > 



We may easily check that C (R) may be defined as a continuous increasing function 
of R £ [0, oo ) such that C(0) = and C(R) — > oo as R — >• oo. Hence, the left-hand 
side of ( 3.17[ ) attains a maximum value for some i?* G (0, oo) and by (3.17) the 
initial bound r can take the largest possible value given by the left-hand side of 
(3.17) for R = R*. In particular, (3.17) may be viewed as a restriction on the 
initial bound r which amounts a restriction on || po II oo and ||mo II oo 3 assuming given 
Sq ■ We also verify that the initial bound can be taken as large as we wish provided 
that 7 — 1 is sufficiently small. 

Now we proceed to prove the compactness of the sequence of approximate solu- 
tions U h . The proof is based on the general analysis carried out by DiPerna in [4] 
and we are going to apply the compactness result in and its extensions in pQ, 
[S] and [7], which together cover the whole range 7 > 1. 

Now, let V h = {p h ,m h ) and F h = (m h ,p h (u h ) 2 +p{p h ,S h )). For any G 
C5°(IR 2 ) we have 
(3.18) 



V h 4> t + F\ 



x[0,oo) 



V h 4> t + F'\ 



■ dx dt 



lV h (jh)Mx,jh) dx 
lV h {jh)\<t>{x,jh)dx- I V h {x,Q)ct>{x,Q)dx 



where 



■ dx dt — 

3=0 

OO 

= E 

3=0 

= E 

3=1 



[V h (jh)j = V h (x,jh - 0) - V h (x,jh + 0). 
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Further, if (rj, q) is an arbitrary entropy pair for ( 1.1 )— ( 1.2 ), with S constant, we 

have 

(3.19) 

oo r (j + l)h 



v h 4>t + q h 0x dxdt^y^ / ^Vt + q h (f>x dx dt 

X[0,oo) j=0^ h J R 



h (x,0)cl)(x,0)dx + V / {q h {jh)lcj){x 1 jh)dx+ / S{(f))dt+ / C{cj>)dt, 

where, for reasons of brevity, we write ry ft = '^(V' 1 , S h ) and g' 1 = <J[(V , S ). Here 

foWl = - 0) - r, h {x, jh + 0), 

and <S(</>) is defined as 

shocks 

[V 1 ] = V h (x(t)-0,t)-r ] h (x(t)+0,t), 

where the sum is over all shock discontinuities (x(t), t) at time t, s = x'(£) denoting 
the shock speed, while C{<f) is defined as 



contact 
discontinuities 



with sum running over all contact discontinuities (x(t),t) at time t, where u h is 
the velocity. The latter is defined over a vacuum interval as the arithmetic mean 
between the velocity at the end of the 1-rarefaction wave bounding the vacuum 
interval on the left-hand side and the velocity at the beginning of the 2-rarefaction 
wave bounding the vacuum interval on the right-hand side. 

We recall that if (77, q) is a convex entropy pair for the isentropic system ( 1.1 1 



(1.2) where S is constant, then 

(3.20) afofc] - lq h ] > 0, 

across each shock wave. Since S h is constant across waves of the first and third 
family, inequality ( 3.20 ) also holds here. Therefore, for any weak entropy pair (77, q), 
we find that the functional 

poo 

S((f>) dt 

is a (signed) measure with locally finite total variation, as a consequence of Re- 
mark o 

Concerning the functional 



C{4>) dt, 

if (77, q) is a smooth entropy pair, we have, in view of previous calculations, 
\u h h h {jh)\ - {q h Uh)j\ < C v e-* h h, 

and so 



< C v diam(K)\[ 



where K is any compact containing the support of </>, which gives that this functional 
is also a measure with locally finite total variation. 

Observe that the weak entropies may be also written as 

r,(p, u)=p ^ <7 (™ + ze^ s ' 2 V 7 " 1)/2 ) (1 - z 2 )\ dz, 



L°° SOLUTIONS FOR POLYTROPIC GAS WITH DIFFUSIVE ENTROPY 



IS 



while a similar formula holds for q. In particular, r\, q are Lipschitz up to vacuum 
if g is smooth. 

We also observe that for the special entropy pair (77*, q*) we have J °° C(<f>) dt = 0. 
Also, for this entropy pair, for nonnegative <j> € Cg^M 2 ) we have 
(3.21) 



E / lvHjh)U(xJh)dx 

00 „ 

= EE / (^ h (^-0)),S^,j/i + 0)) 

- 77* (^(z, j/i + 0), jTi + 0))) 4>{x, jh) dx 

00 „ 

"EE / (v*(y h (xjh-o)),s h (x,jh + o)) 

j=i iez Jli 

- n * (v h {x,jh - o), s h {x, jh - o))) 4>{x, jh) dx. 



The first sum in the right-hand side of equation (3.21) is nonnegative for nonnega- 



tive <j>, since V h (x,jh + 0) is the average of V h (x,jh — 0), in each interval Ii, and 
77* is convex. Therefore, we get 



(3.22) 



E / H(jh)U(x,jh)dx 

00 . 

>-EE/ i t V-)(s l (i 1 i'> + o)-s k (^j/ 1 -o)Mi ) i# 

00 . 

> - E Ce^' l /i / 0(x, jh) dx, 



where ■ • ) = Jq V*s(V h ( x 'jh — 0), .4(0)) is the coefficient of the linear re- 
maining term in the trivial Taylor expansion of zero order in the variable S and 
A{6) = (l-9)S h (x,jh-0)+9S h (x,jh+0). In particular, both the left-hand side as 



well as the second term of the right-hand side of ( 3.21 1 are measures of locally finite 



total variation. As a consequence, we may apply equality (3.211 with cj> replaced 
by the characteristic function of any suitably chosen rectangle x\ < L = MAx, 
< t < T = Nh, to find that 

(3.23) E E / D 2 v r ] ':(---)(V h (x,jh-0))-V h (x,jh + 0)) 2 dx< const., 

jh<N \iAx\<M h 

for any M, N > 0, the constant depending on M, N, where D v r]%(- ■ ■ ) = J (1 — 
9)D v r)*(B(9), S h (x,jh + 0)) dd is the coefficient of the quadratic remaining term in 
the Taylor expansion of first order and B(9) = (l-9)V h (x,jh+0)+9V h (x,jh-0). 

Since for all weak entropy 77 we have l-Dy^l < C^Dy-q^, for some C v > 0, it 
follows from (3.231 that 



(3.24) 



E E / \DvV\(V h (x,jh-Q))-V h (x,jh + Q)) 2 dx 

AT 1 „- \ ~ I S l\ r •> h 



jh<N \i&x\<M' 

for any M, N > 0, the constant depending on M, N. 



< const., 
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We can then use DiPerna's method in [3] to prove the W lo *' 2 compactness of the 
distributions rj' t l + q% by decomposing the functional 



m = Y] / lv h (jh)Mx,jh)dx 

3=1 Jit 



as 



(3.25) 



oo - 

iW = EE / lv h Uh)Mx,jh)dx 
3=1 i& J Ii 

oo - 

= EE(^'^) / lv h (jh)jdx 

.7 = 1 l& Jl * 

lr, h {jh)\(4>{x,jh) - <t>{x h jh))dx 



Li(0) + L 2 (0). 



We consider the two terms separately. We have 



oo „ 

£iW>)=EE^'^) / lv h {3h)lv + h h {jh)\sdx=:L ll { ( t>)+L l2 {4>), 
7=i iez 



where, if [ij(V, S)\ = r)(V-,S-) - r)(V+,S+), we denote 



Since U^O'/Olsl — ^ e we dearly have 



\LM<f>)\ < cu\\ 



Concerning Ln(<j)), we have, cf. (3.24), 
(3.26) 



\LiM\ < 



3=1 iGZ 



oo - 

EE^ 1 '^) / Dfa h (---)(V h (x,jh-0)-V h (x,jh + 0)) 2 dx 
7=1 iez 



3=1 i6 

< C 1101 



Hence, we have 



|£i(0)l<CiM 
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Next, exactly as in 0], we find, assuming that the test function <f> satisfies 
supp0 C [-N, N] x [—J, J] and keeping 9 > sufficiently small, 
(3.27) 



M0)i< EE/ HWiw^jh)- 

\j\<J\i\<N Jli 

E / I^WllA^dx 



<f>(i,jh)\ dx 



< 



< 



< 



< 



< 



\3\<J 
\i\<N 



\j\<J 

\i\<N 

i-E E 

|j|<J|i|<JV 

(Ax 2a+1 



Ax e 



ix e [ \lv h (jh)]\ 
Jii 



dx 



■lie- ( 



Ax e 



(2J + l)(2iV+l) + A a ; e £ / |[V l (^)l| 2 dx) 
b'l<J R 



\i\<J 



1011c ( 



Aa; 2a+1 



A 



-2 N^Nca 

where C Q denotes the Holder space with seminorm 



= sup \<f>(x) - </)(y)\ I \x - y\ a , a > 1/2, 



and where C 2 depends on the support of <fi. Thus 

IM0)l < Cill^lloo, and \L 2 (^)\ < C 2 (Axf\\<t>\\c« 

for appropriate a,j3 £ (0,1), for some positive constants C\,C% depending on 
supp <j), but independent of <f>, and through the Sobolev imbedding theorem 

Luffl < C^AxfUWwx,,, 

for an appropriate q € (1, 2) and constant depending on the support of (p. 

In this way we obtain by the usual interpolation argument that for any weak 
entropy pair (r],q) for (1.1 )— ( 1.2 ) we have 

V (V h , S h ) t + q(V'\S h ) x £ { compact of W^' 2 (R x [0, oo)) }. 

We can then use the compactness results in [TJ [H [7] to deduce that we may 
extract a subsequence of (p h ,m h ,S h ) converging in _L 1 1 oc (]R x [0, oo)) to a weak 
solution (p(x, t),m(x, i), S(x,t)) to ( 1.1 h |L5| ). Also, (3.22) implies the entropy 
inequality (1.15), and (3.18) implies (1.10) by a calculation similar to the estimate 
for L2(4>) above. 



Concerning the decay property (1.16), we prove it as follows. First, from the 
above discussion, we deduce that for any weak entropy pair we have 

\(il(p,m,S) t + q(p,m,S) x ,(l>)\ < CiH^U, 

with C\ depending only on suppc/) and bounds for (p,m,S). Hence, if U T = 
(p T ,m T , S T ) is the self-scaling sequence U T (x,t) = U(Tx,Tt), we see that for any 
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entropy pair 

\{rj(p T , m T , S T ) t + q{p T , m T , S T ) X , <f)\ < Ci||0||oo, 
while from we have, for < t < T, 

~/*(p, m, S)(x,t) dx > / rj*(p,m 7 S)(x,T) dx — C I I e~ s dxds 

[0,L] J[0,L] Jt J[0,L] 



> / ?~*(p, m, S)(x,T) dx — CLe * 

J[0,L] 



Hence, we can apply the decay result in [2] to deduce (1.16), which then concludes 
the proof. 
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